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I. INTRODUCTION 


In recent years laser radar (lidar) systems have developed Into useful 

tools for remote monitoring of the earth’s atmosphere. Satellite-borne 

lidar systems provide the attractive feature of wide coverage and in some 

cases may be the only practical means of probing the upper atmosphere. 

The quality of lidar data depends of course on the system noise characteristic 

The receiver output will be contaminated by the usual background and shot 

noise components and in some cases by speckle noise. In certain situatioi s 

speckle noise dominates and can seriously limit system performance. 

Laser speckle has been studied extensively. Goodman* published one 

of the earliest papers describing speckle effects on optical radar 

performance and recently an entire issue of the Journal of the Optical 

2 

Society was devoted to speckle. In this report we develop and summarize 

the equations which relate the statistics of the speckle noise at the 

receiver output to the lidar system parameters. Fortunately, we were 

able to adapt much of the existing theory to the satellite-based lidar 

problem. Much of the material was obtained from the books Laser Speckle 

3 

and Related Phenomena edited by J. C. Dainty and Statistical Properties of 

4 

Scat tered Light by Crosignani, Di Porto and Bertolotti and from the 
November 1976 special issue of JOSA on speckle. 

The lidar system model is described in Section II, and in Section III 
the statistics of the signal and noise at the receiver output are derived. 
Scattering media effects are discussed in Sections IV and V. Polarization 
and atmospheric turbulence are considered in Sections VI and VII. And 
finally, in Section VIII, the major equations are summarized and evaluated 
for some typical system parameters. 




II. SYSTEM MODEL 


A typical satellite-based lidar system is illustrated in Figure 1. The 
transmitting telescope is used to project an image of the optical source 
onto the scattering medium. In most cases the source would be a narrowband 
laser. However, since speckle noise is significantly influenced by source 
temporal coherence, we will assume only that the source coherence bandwidth 
is a small percentage of the center wavelength. Although a major portion of 
the propagation path will be in free space, in the lower atmosphere turbu- 
lence effects could become important and will be considered. The scattering 
medium could be a rough surface such as the earth or an ensemble of 
scattering centers such as aerosols (clouds) and air molecules. 

The receiving telescope model is illustrated in Figure 2. The spatial 
filter following the objective lens is normally a f ield-of-view (FOV) stop 
which is used to reject background radiation. After the field stop the 
optical signal is collimated, passed through an interference filter and then 
focused onto a pnotodetector . Finally, the signal from the photodetector 
is passed through an electrical filter to further limit the system 
bandwidth. 

The problem we are concerned with is the detectability of the scattered 
signal. The receiver output will be contaminated by the usual background and 
shot noise components and in some cases by speckle noise. In certain 
situations speckle noise dominates and can seriously limit system performance. 
In the following sections we will present the equations which relate the 
statistics of the photodetector output to the lidar system parameters and the 
characteristics of the scattering medium. 
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Fig. I. Lidar Geometry 


/ 



Objective 


Fig. 2 . Receiving System Model 





III. STATISTICS OF THE RECEIVED SIGNAL 


The signal statistics at the receiver output can be deduced from the 
statistics of the optical radiation incident on the receiver aperture. Let 
u (_r , z , t ) be the analytic signal representation for a single polarization 
component of the electrical field. For a quasi-monochromatic field u takes 
the form 

u(_r,z,t) = |A(£,z,t) |e iQ e iu) c t (1) 

where 0 is the phase and u> c is the center frequency. For the moment we will 
assume the scattering medium is diffuse so that the received field consists 
of a large number of randomly phased contributions arising from different 
scattering centers. *f the scatterers are randomly distributed over distances 
which are large compared to the optical wavelength, the scattered signal will 
be incoherent. Under these conditions A is Gaussian distributed while the 
intensity obeys negative exponential statistics^ 

I(r,z,t) - | A| 2 (2) 

P (I) = 7F ex P<- TTT* (3) 

where <I> is the average intensity. 

The scattering conditions can be relaxed to include the cases where the 
surface is not rough compared to A and where the scattered field contains a 
coherent specular component. These situations are discussed in Section IV. 

To calculate the signal statistics at the receiver output.it will be 
necessary to evaluate the mean and covariance function of the intensity 
flue tuations 
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C I ( -l» t l’-2’ t 2’ z) = <I (l 1 » z » t 1 ) I (£ 2 ,z,t 2 ) > “ <1 (£ i » z » t i) > < I(r 2 ,z,t } )> . 

(4) 

Since the field A is a circular complex Gaussian random variable, the 
covariance function can be expressed in terms of the mutual intensity of 
the field^ 

C X = | <A(r 1 ,z,t 1 )A (£ 2 ,z,t 2 )>| t = I J A (-i» t i»- 2 » t 2 » z ) I * * 

The mutual intensity function depends on the characteristics of the 
scattering medium and the incident optical signal. In Sections IV and V 
we show that for most casts of interest it can be assumed that A is spatially 
homogeneous in the transverse direction. The temporal dependence of A arises 
from a number of factors. Movement of the scattering centers with respect 
to the receiving telescope caused by Brownian motion, turbulence, wind, or 
spacecraft motion will introduce temporal fluctuations. Those effects are 
generally stationary. However, if the optical source is pulsed and the 
density of scatterers varies with distance (z) , the average value of the 
received intensity will change with time. This effect can be included in 
our analysis by writing the mutual intensity function in the following form 
(see Section V) 

J A^-l' t l*-2* t 2 ,Z ^ = D ( t i* t 2 ) ^A ( -l " -2^1 _ 
r>(t ,t 2 ) is given by 

CO 

D( t lf t 2 ) - / dzP(t 1 - ~)P(t 2 - ~) -f" (7) 

where P is the incident pulse shape and p is the number density of scatterers. 
If the scattering medium is a rough surface, p Is a Dirac delta function and 
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s 


becomes 

3 2 

Dftj.tj) - P(t ) P(t 2 -^) -§ (8) 

z 

where (3 f is the surface reflectivity. Using Equations (6) and (7), the 
intensity covariance function can be written as 

C I = D2(t i’ t 2 ) " -2* t l “ t 2 ) | 2 * (9) 

The average intensity is obtained by simply evaluating the mutual intensity 
function at the point r^ = jr^ and t^ - 

<I> = D(t,t)J A (0,0) . (10) 

We assume the receiver FOV and interference filter are adjusted to 
admit all the scattered signal energy which is incident on the telescope 
objective. The total received signal intensity can therefore be written as 

I g (t) = /d 2 rW(r - vt)I(r,z,t)/ /d 2 rW(r) (11) 

where W is an appropriate real and positive aperture weighting function and 
v is the transverse spacecraft velocity. W is dimensionless. Note for a 
point detector W is just a Dirac delta function. 

Two types of receiver structures will be considered. The analog 
receiver is applicable to strong signal conditions where h(t) is a low-pass 
filter (see Figure 2). When the signal is weak, h is an integrator whose 
output is proportional to the number of detected photons. In both cases the 
noise at the receiver output will depend on the statistics of 1^. 





The mean and variance of the signal at the receiver output can be 
calculated using Campbell' 3 theorem^ 

E[s(t) ] - o<I g (t)> *h(t) (12) 

Var[s(t)J = a <1 (t)> *h 2 (t) 
s 

ao no 

+ a 2 / dtj^ / dx 2 C I (T lt T 2 )h(t - Tj)h(t - t 2 ) (13) 

— ao — ro S 

•./here a is a constant which depends on the photodetector efficiency and 

c 

is the temporal covariance function of the received signal. 

The variance of S can be regarded as the signal induced noise. The 
first term in (13) is the quantum or shot noise component arising from the 
random emission of photoelectrons by the detector. The second term arises 
from the speckle induced fluctuations in i s< 

The average value of I is calculated from (11) using (10) 

<1 (t)> = D(t,t)3f.(0,0) . (14) 

S A 

The temporal covariance of I can be written in terms of J. and D 

S A 

C I = D ^ t l ,t 2^ I d2 ll^ A (l+ VT,T) | 2 r w (c) (15) 

s 


T 


- t. 


R is the autocorrelation function associated with the receiver 
w 

aperture weighting function 

(jr ) = /d 2 rW (£)W*(p + r) . (16) 

The average value of S(t) can now be evaluated by substituting (14) 
into (12) and using (7) 

00 

E[S(t) ] = a J. (0,0, z) / dzP*" (t - — )*h(t) . (17) 

A 0 C 
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This expression can be further simplified by assuming the transmitter pulse 

width and receiver bandwidth are chosen so that the small scale structure 

2 

in p(z) can be resolved. In this case p(z)/z is approximately constant 
over the important range of integration in (17). By choosing the time 
reference so that P and h peak at t ■ 0 we have 


E[S(t) J = oj (0,0,2) £ - ( — ^ J dzP 2 ( t - ^r)*h(t) . 

A (ct/2 r o c 

Under the same conditions the expression for the variance becomes 


(18) 


where 


Var [s(t) ] = aJ.(0,0,z) j dzP 2 (t - — ) *h 2 (t) 

A (ct/2) Z 0 c 

+ « 2 t I* 2 ?- i® dT lV- + - T * T) l 2R p (T)R h (T)R w ( i ) 


R (x) = / dtP(t)P(t + t) 
P i® 


(19) 


( t ) = / dxh(t)h(t + t) 


Rp and R^ are the autocorrelation functions associated with the transmitted 

pulse shape and impulse response of the receiver filter. 

It is instructive to consider the physical implications of (19). If 

the receiver employs a point detector of infinite bandwidth, both R^ and 

8 are Dirac delta functions, in which case the second term in (19) reduces 

to the square of the average intensity. For a finite aperture and finite 

bandwidth receiver, R and R. are rather broad functions of r and t with 

w h — 

'V, 

peak values at r = t = 0. Since J. and R are also maximum when their 
- A p 

arguments are zero, (19) will be less than the average intensity when a 
finite aperture and finite bandwidth receiver is used. In this case the 
receiver spatially and temporally averages the intensity fluctuations. The 
averaging will depend ipon the relationship between the spatial and temporal 
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coherence lengths of the received signal and the receiver aperture size 

and bandwidth. The averaging will also depend upon the transmitted pulse 

2 

length through the factor in (19). At any given time t, the signal 
energy incident on the receiver aperture arises from an Illuminated 
scattering volume whoLe thickness is determined by the transmitter pulse 
length. The scattered signals arising from each of the illuminated 
scattering centers add on an amplitude basis. Consequently, the spatial 
statistics of the intensity are unaffected by the pulse length (aside from 
a possible scaling factor). However, since the signals scattered by different 
centers are uncorrelated, the received intensities at two different times 
will be purtiall/ correlated only if the corresponding scattering volumes 
have some points in common. Overlapping will occur only if the time 
difference is smaller than the transmitted pulse length. The intensities 
will be uncorrelated if the time difference is larger than the pulse length. 


The temporal statistics are also affected by the source coherence, random 
motion of the scattering centers and spacecraft motion. These effects 
appear in the .If factor in (19). The spatial statistics are determined by 
the '^ansmitter beam divergtnee, characteristics of the scattering medium 
and propagation t 'fects such as turbulence. These factors are discussed in 
more detail in Sections IV and V. 

The relative importance of spatial and temporal coherence length, pulse 
length, aperture size and receiving bandwidth can be determined by evaluating 
(18) and (19) for the idealized case where all the functions are Gaussian 



|J A (0.0) 


exp (- 


2e 2 

c 



( 20 ) 


P ( t ) 



( 21 ) 
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h(t) - exp(-t 2 B 2 ) 


( 22 ) 


2 

f 

W(r) ■ exp( 

R 


where 

p = transverse spatial coherence length of signal 

c 

t * temporal coherence length of signal 
c 

t « transmitter pu’se length 
B » detector bandwidth 
R *> receiver aperture radius. 

Using (20) through (23) in (17), (18) and (19) gives 


E[S(t)l 


o3f A (0,0) 


P (ct/2) 
(ct/2) 2 



Var|S(t)] 


Mi) . MS ). 2 

Ji 


(23) 


(24) 


(25) 


where 


M„ - J + 


1/2 


«T- 


1 + 


1 

B 2 , 2 

c 


bV 

p 


2 2 2 
B*(R* + p*) 


(26) 


(27) 


The equations are valid for scattering by tenuous media such as low 
density aerosol distributions and air molecules. For rough surface 
scattering the 

p (ct/2) 

(ct/2) 2 


factor must be replaced by 



2 

7 * 
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ii (24) and t set equal to infinity in the expression for M . 

P 1 

The average signal varies with time in proportion to the density of 

_2 

scatterers. The (ct/2) “ factor accounts for decay in signal strength as 
the distance between the spacecraft and illuminated scattering region 
increases with time. The first term in the variance expression is the 
shot noise component. It is proportional to the mean signal. The speckle 


noise component is proportional to the square of the mean signal and 
inversely proportional to MgM^,, the effective number of received cor- 
relation cells. Mg is the effective number of spatial correlation cells 
averaged by the receiver aperture while is the effective number of 
temporal correlation cells. 

To reduce speckle noise we would like to make Mg and ftj, as large as 
possible. The parameters R, B and V are generally inflexible. The 
aperture size R is constrained by weight limitations and the physical 
construction of t he spacecraft. The detector bandwidth B is determined by 
the required range resolution. And of course the spacecraft velocity is a 
function of the orbit. Thus the only parameters which may be adjusted are 

p , t and t . 

c c p 

The coherence time of the scattered radiation depends on the scattering 

medium and source coherence. If the scattering medium is a rough surface 

then t is exactly equal to the source coherence time. If the scattering 

mechanism is resonant absorption and reradiation, then r is equal to the 

inverse Doppler line width. For air molecules dominated by Brownian motion, 

i is on lie order of 1 ns. If Rayleigh scattering from air molecules or 

-2 -2 - 1/2 

aerosol and particulate scattering is involved, then i = (t , + t ) 

c d c 

where is the source coherence time and 1/tj is the Doppler line width, 
tj is on the order of 1 ns for air molecules governed bv Brownian motion and 
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1 ms for aerosols or particles embedded in a turbulent flow (e.g., clouds 
and smoke plumes). Therefore, except for resonant scattering, can be 
controlled by varying the source coherence time. 

Mj, can also be increased by reducing the source pulse length . But, 
there is a practical limit to the reduction. To keep the per pulse signal 
energy constant, the peak intensity must be increased as i is decreased. 


However, when resonant scattering is involved, the peak energy density in 
the scattering medium must be maintained below saturation levels. This 
problem is most important v/hen the density of scatterers is low. 

In the absence of strong turbulence effects is a function of the 
size of the illuminated spot on the scattering medium and the distance 
from the scatterers to the receiver. If the whole transmitted beam illumi- 
nates the s . tering medium we have (see Section IV) 


where \ is the wavelength and 0 
the transmitter is operating at 


P c = 0: 


(28) 


is the transmitter divergence angle. If 


the diffraction limited divergence, then 


9 


T 



(29) 


where R^. is the transmitter aperture radius. Usually the transmitter and 
receiver apertures are nearly equal in size. Thus Mg = 2 for diffraction 
limited operation. The transmitter divergence normally will be much larger 
than the diffraction limited value to minimize pointing and tracking problems. 
If the transmitter divergence is increased, the receiver FOV must also he 
increased so that the receiving telescope sees all of the illuminated 
scattering volume. Unfortunately, background noise power increases witn 
FOV. If the background noise results from an extended diffuse source, its 
characteristics will be <£.nilar to the scattered signal since both are 
incoherent and Gaussian. 
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Uniform background noise sources can be described by their spectral 
radiance function, N(f), which is defined as the power radiated at 
frequency f per hertz bandwidth into a unit solid angle per unit of source 
area. 7 If the source fills the receiver FOV, the effective average noise 


intensity incident on the receiver aperture is 

<I n (r,z,t)> = J n (0,0) = N(f)B r n r . (30) 

is the noise mutual intensity function, is the bandwidth of the 

receiver's interference filter and 0 is the solid angle FOV of the receiver. 

r b 

Because fl is matched to the transmitter divergence and raav change, it is 
convenient to normalize I with respect to a standard FOV such as the 
diffraction limited value. 


R0 

<I n > = ■'nd l (0 - 0)( - T ' ) 


J nd« (0 - 0) = N(f)BjO^) 


(31) 


J is just the effective noise intensity for a diffraction limited FOV. 
ndi. J 

The mean and variance of the noise voltage at the receiver output can 
easily be calculated if we use the form in (20) for the noise intensity 
correlation. 


R0 t _ 1/2 Re ‘ 

Eln(t)] » E(n dJ )(^I) - -J^O.O) - (-^J 


(32) 


Kn dl > 2 . «"„> 2 «• * 

Var|n(t)] + — — - , 


(33) 


✓2 


ST 


where 


R0 t 

N S = i + <— > 


( 34 ) 
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1/2 


(35) 


N„ = 


B 


1 + 


9 2 2 2 2 

b (r + x z /ep 


The spatial coherence length of the background noise is determined by the 
receiver FOV and is equal to \/Q^. The temporal coherence length is 
determined by the interference filter bandwidth and is equal to 1/B { . 

Background noise contaminates the receiver output by adding a DC com- 
ponent and introducing additional shot noise and speckle noise. Normally 
the speckle noise component will be small since the interference filter 
bandwidth is usually much larger than the detector bandwidth. For example, 
the bandwidth of a 1 % interference filter centered at 530 nm is nearly 
100 CHz compared to detector bandwidths of 1 GHz or less. Thus N^. is on 
the order of 10^ or larger. 

Although signal induced speckle noise will decrease with increasing 
0^, background shot noise and speckle noise will increase. The optimum 
transmitter divergence can be determined by minimizing the sum Var(S) + 
Var(n) with respect to 0^. 


1/4 


optimum 


E(S/ /M t 

E(n d( )/^ + E(n dl ) 2 /N T 


(36) 


If the optimum transmitter divergence is used, the total noise power is 

1/2 


„ „ E(S) . 2E(S) 

Var(S) + Var(n) = — » — + --rrr- 

JZ .. 1/2 


E(n dt ) E<n di ) ' 


/2 


N, 


. i + 

Jo * 

■ 2 L 


2 5/2 K(n d( ) 

1/2^ 

M 

T 


L 1 



(37) 
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E(S) and ECn^) may be viewed as the average number of signal and noise 
photons detected in a time interval of length 1/B. Usually Efn^) will be 
very small since it applies to the diffraction limited receiver. In most 

cases, therefore, the total noise in (37) is approximately equal to the 
signal shot noise. However, it may not be practical to achieve the opti- 
mum value in (37), particularly if the transmitter divergence given by 
(36) exceeds the angular size of the target. 

Although the mean and variance of the signal plus background noise 
are useful parameters for accessing the lidar system's performance, a 
complete analysis requires knowledge of the signal and noise probability 
distributions. The instantaneous intensity evaluated at a single point is 
exponentially distributed if the scattered signal and noise are both 

completely incoherent [see Equation (3)] and the photon counting statistics 

4 

are Bose-Einstein distributed. Because the receiver integrates the 
intensity over space and time, the intensity statistics will be modified. 
Goodman^ has shown that the gamma distribution is a good approximation for 
the probability density of the integrated speckle patterns 


P ( I) = 


<7?~) i"" le *p(- m 717) 

r (m) 


(38) 


I > 0 


where m is equal to M^M^, for the scattered signal and N^N.^. for the background 
noise. The photon counting distribution associated with (38) is the 
negative binomial distribution. If desired, the exact probability density of 
integrated speckle can be calculated by expanding the scattered field in a 
Karhunen-Loevre series.^ 

The signal statistics derived in tMs section are summarized in 
Section VIII where they are evaluated fot some typical system parameters. 
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IV. SCATTERING BY A ROUGH SURFACE 


In the previous section it was shown that the speckle noise power 
could be calculated from knowledge of the intensity covariance function 
of the received signal. 

Cj « <I(r 1 ,z,tj) I(r 2 ,z,t 2 )> - <I(r L ,z,t 1 )><I(r 2 ,z,t 2 )> . (39) 

If the surface is rough compared to the optical wavelength, the scattered 
field is a complex circular Gaussian process. The magnitude of the 
scattered field is Rayleigh distributed and the intensity is exponentially 
distributed. In this case the intensity correlation function can be 
written in terms of the mutual intensity of the field."* 

Cj = |<A(r ]t z,t 1 ) A*(£ 2$ z,t 2 )>| I 2 = | J A (£ 1 ,t 1 ,£ 2 ,t 2 ,z) | 2 . (40) 


If the surface also possesses a specular scattering region which is 
too small to be resolved by the receiving telescope, then the scattered 
field will contain a coherent component. In this case the intensity 
distribution is a modified Rician density."* 


P(I) = 


<1 > 

r 


exp 


I + I 


L Av 


<1 > 

r 


! 0 2 <1 > 


I > 0 


(41) 


where 1^ is the random intensity component and I ^ is a coherent specular 
component. is a modified Bessel function of the first kind. When the 

rms fluctuations of the surface roughness are not large compared to the 
wavelength, the surface can be decomposed into two identically shaped 
surfaces which are colocated in space, one of which is very rough and the 
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other of which is very smooth but with a reduced reflectivity. In this 
case the intensity distribution is also a modified Rician density. 

The mutual intensity in the observation plane can be written in terms 
of the mutual intensity in the scattering plane using the Huygens-Fresnel 
principle. If the observation plane is in the Fraunhofer zone, we have 


A(r , z, t) 


Xi exp 



d pA (p,t) exp 
- s ~ 



(42) 


where A^ is the scattered field directly in front of the scattering plane. 
In most cases the scattered field can be approximated by the formula^ 


A H (p,t) P r <V(p - vt,t) exp 



(43) 


6 r is the average surface reflectivity, A^ is the complex amplitude of the 
incident field, v is the transverse velocity of the spacecraft and h(p) is 
the random height distribution of the scattering surface. Equation (43) 
is essentially the geometric optics approximation for the scattered field. 

It is reasonably accurate provided the surface slopes are small. The 
formula implies that the scattered field is reduced in amplitude and phase 
shifted with respect to the incident field. 

The mutual intensity of the scattered field is calculated by assuming 
the surface height fluctuations are a zero mean homogeneous Gaussian process 


J A ( £i-£2 


t 2 > m e r J A ^-1 ~ ~ vt 9 , t Jf t 2 ) 

i 2 

[-(t) t c h (0) - c h<£i - e 2 > 


• exp 


(44) 


where is the mutual intensity of the incident field and C^(p) I s the 
surface height correlation function. The mutual intensity in the observation 
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plane can now be calculated using (42) and (44) 


**a(-1*-2* *■!' ^ 9 ’ ^ 



- t 2* t r t 2 ) 


( 45 ) 


Note that in (45) we have neglected some multiplicative phase factors. 

These factors do not affect the intensity correlation function (Eq. (40)] 
because it depends only on the magnitude of J . 

t\ 

The major contribution to the integral in (45) occurs for near 
If the receiving telescope is unable to resolve an object in the scattering 
plane whose scale size is on the order of the spatial coherence length of 
the scattered field, may be approximated by^ 


r f 2 

= ~ 2 L 2 d £ J a " - (t l " t 2>» t i»t 2 ] 

\ z -L i _ 


• exp 


. 2tt 

1 Tz e • <£i - -2 


(46) 


Now, the incident field amplitude is related to the transmitter 


aperture distribution A^. 


t l . IT 2\f ,2 . , . /. 2 tt \ 

^(p.t) = — exp i-i — p J j d uA^ (u,t) exp I i — u • p j 


(47) 


it is convenient at this point to express A^. explicitly in terms of its 
temporal and spatial dependence. 


Vs-') - 4 - t) T (' - t) 


(48) 
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P is the transmitted pulse shape, T is a random time function which is 
related to the source temporal coherence, and A^, is the spatial field 
distribution in the transmitter aperture. Using (48) and (47) we have 


J A 1 ( ^1^2 ,t l’ t 2 ) 


( C 1 - T) P (h - T) C T (t l - £ 2> 

P [' 1 Tz (*! * "zj (I 

(j If <Hi • fii - h 2 • £2] 


.2 2 
A Z 


d ^ d 2 u 2 A^u^) A t (u 2 ) 


(49) 


where C^, is the source temporal coherence function. Substituting (49) into 
(46) and integrating over p gives 


J A^-l*-2’ t l’ t 2 ,Z ^ = D ( t i* t: 2^ ‘V-l “ -2 ,t: l ” l 2^ 


D( 


6 2 P ft. - — ^ 

t ) = - r - A J ll 

1 z 2 ; 


p it 2 

_AJ £ 


(50) 


(51) 


M T) f 

J A ( - ,T) = 4 2 

A z 


d u A t (u) A t (u - r - vt) exp 


. 2tt 

-i — u • vt 
Az - 


(52) 


We will evaluate (52) by assuming the aperture illumination is Gaussian. 
This is a reasonable assumption if the laser is operating in the funda- 
mental Gaussian mode and if the transmitting telescope does not obscure the 
central region of the beam. The integral in (52) can be simplified by 
noting that the magnitude of u can never be greater than the transmitter 
aperture radius R^. If R^vT/Az << 1, a condition which is usually met, the 
exponential in (52) can be neglected. Therefore, if we let 

f 2 2*1 


A t (u) « A t exp 



(51) 
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J. becomes 
A 


J A ( r,T) 


1 i c t( t ) exp 


r | .2 

/ 2 2 \ 


| r + vt | 

1 . tt R ) 



1 T 


2 

,2 ,2 .2 


L 

x f T y 

md 


(54) 


where is the effective transmitter focal length and Ij is the average 

incident intensity. It is convenient to write (54) in terms of the 
effective transmitter divergence angle. 


.I A (r, T ) 
A - 


1 1 C,j, (t ) exp 


L 2 x 2 J 


If the transmitter 
limited value A/R^ 


6. * Hr 



beam is collimated, f r 


The spatial coherence 


1/2 


= 00 and 
length 


0^ becomes 
is defined 


(55) 

(5f>) 


the diffraction 
as 


P 


c 


A/0 t 


(57) 


Notice that the spacecraft velocity effect is doubled because both 
the transmitter and receiver are moving. The variance of the received 
signal is evaluated by integrating J^(r 4- vT,T) over the receiver aperture 
[see Eqs. (14) and (16) J. The vT parameter was introduced bv the receiver 
motion. Transmitter motion Introduces the vT parameter in Eq. (55). Thus, 

v should be replace by 2v in the analyses in Section III. 

% 
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V. SCATTERING BY A RANDOM DISTRIBUTION OF SMALL PARTICLES 


Scattering by an ensemble of particles differs from rough surface 

scattering because the particles are generally in motion, resulting in a 

continually changing scattering medium. The trajectories of molecules 

are usually governed by Brownian motion while the motion of material 

particles such as aerosols are determined by the surrounding medium. 

The book, Statistical Properties of Scattered Light , by Crosignani, 

4 

Di Porto and Bertolotti, gives an excellent treatment of this subject. 

ihe spatial statistics of the scattered field are very similar to 
those for rough surface scattering. However, the temporal statistics 
are considerably different because of the random particle motion. The 
evaluation of the scattered field is greatly simplified by using the 
Ray leigh-Gans hypothesis. This hypothesis is similar to the weak 
scattering assumption often used in turbulence problems. Physically, 
the Rayleigh-Gans hypothesis implies that the scatterer introduces a 
small optical perturbation to the surrounding medium and that each point 
of the scatterer sees the incident field practically unperturbed. Under 
these conditions the field scattered by the n t * 1 particle can be written 
in tne form 


A * A, (p - vt,t,z) R (t) exp [ ik • r ( t ) ] . (58) 

sn l -n — n - -n 

Aj is the complex amplitude of the incident field, (p^.z) are the 
particle coordinates and k is the optical wavenumber vector, r is the 
trajectory of a given point on the particle and R ^ depends on the rota- 
tional motion around r . When the particle is spherical and r is the 

-n -n 

trajectory of the center of the sphere, R is constant. 

If the observation plane is in the Fraunhofer zone, the field 

L 

contribution from the n scatterer is given by 
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(59) 


A n (r,z,t) “ 


sn 

Az 


exp 


( . it 2 2tt \ 

[- 1 Ti r 1 7z ‘ -) 


The total field in the observation plane is calculated by integrating 
(summing) (59) over all the illuminated particles. Since particles may 
also be distributed along the direction of propagation (z), the integral 
must extend over both the longitudinal and transverse dimensions. 


A(r,t) 


r 

dz 

J 0 


sn 


d 2 p 

-n Az 


exp 



i 


2tt 

Az -n 



(60) 


Note that the integral notation in (60) is used for convenience. In reality 
the total field is calculated by summing over all the illuminated scatterers. 

Equation (60) can be simplified by noting that at any given instant 
the transmitting pulse illuminates a relatively small region of the z axis. 
Therefore, the phase factor exp -i ~ r 2 is essentially constant over 

A Z 

the important range of integration. It can be pulled outside the integral 
signs and dropped since we are only interested in the magnitude of J^. 


A(r , t) = 


dz 


9 A 
.2 sn 
d p “T — 

-n Az 


exp 


(‘ fi £n ‘ i) 


(61) 


The mutual intensity function for the scattered field in the observation 
plane can now be calculated using (58) and (61) 


r 


‘ 1 A ( -l*-2 ,t l ,t 2 ) 


dz 


dz 


-n 


. -A A* > 
,2 sn sm 
d p 


A z z 

n m 


-m 


* exp 


! X 2* I 2” 

i — P • r - i - — p • r 
> z -n -1 Az -m -2 
n m 


(62) 




where 



<A an A * > " <A ,(P n ~ Vt ,t ,Z ) A *(p " Vt,,t,,Z )> 
an sm 1 -n -lln 1 -m - L i m 


• <R (t.)R*(t 9 ) exp {ik • [r (t.) - r (t„)]}> . (63) 

n i m i - -n l m i. 


Whenever the rotational motion is decoupled from the translational motion 
(this is not always the case), (63) becomes 


<A A* > - <A A*><R (t )R*(tJ><exp { ik • [r (t ) - r (t 7 )]}> . (64) 

so sm i i n l m 2 - -n 1 -m 2 


The ensemble average of th^» exponential term in (64) is taken over the 
initial positions and velocities of the particles and over fluctuating 
quantities of the medium that influence particle motion. Equation (64) 
can be further simplified if we assume the medium is stationary and 
homogeneous, which implies that the initial parricle positions are uniformly 
distributed 


<A A* > * <A . A*>«"R (t. - t n )R*(0)> 
sn sm i i n 1 2 m 


• <exp {ik • fr (t. - t,) - r (0) ] }> . (65) 

- -n i L — m 


The particle trajectory is given by 


*n (t) " r n0 + J 0 dC -n^ , -n0’- n 0 ) 


where 


v(t> " 77 r(t) 
-n dt -n 

JnO * In (0) 

?n0 ' !!n (0> • 


Using (66), the average of the exponential factor in (65) now becomes 
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<exp { ik • (r (t) - r (0)]}> - 
- -n -m 


<exp i l b • <r n0 • in.o )! exp 


ik 


o ^ v n (t,r n0 ,v n0 )j> . 


(67) 


Because of homogeneity the right-hand side of (67) must be independent of 
the initial positions of the n^' and m 1 *' particles. This will occur only 
if the expectation vanishes for n i m. Under these conditions (62) reduces 
to 


(t j » r 2 » f jtt^z) 


fa, f 

0 J 


■ <R n (t l 


d 2 p r Aj (p - vtj.tj) A ( (p - vt 2 .t 2 )> 

A z 

- t.,) R*(0)><exp (ik • (r^tj - t 2 > “ ^nO^ 

(68) 


where p(z) is the number density of scatterers. We have assumed that the 
scattering density is essentially constant throughout the incident beam 
cross section. 

Further evaluation of (67) follows directly f rom the work on rough 
surface scattering presented in the previous section. The last two 
factors in (67) are independent of p and z and can be removed from the 
under integral signs. If in addition we use the representation in (47) 
and (48) for the incident field, (68) becomes 


J A<-r-2*t i*t 2 » z ) “ D(tj.t 2 ) ^ A (r,0 


D(tj ,t 2 ) 


dz P t 


0 


, - •:■) <(>,- :) 


2S) 

o 


3^(r, t ) - C t (t) C (i) Cj.(i) -^2 f cTuAj(u) A.j.(u - r 

* 7 # 

r . 2 , 1 

• exp | -i — u • vt | 


" VT) 


l 


where 


C R ( ° ■ ' R n <T > ^ <0> ' 

C p ( r ) - 'exp ( ik • | r ^ ( t ) - r^ | ) * 


(69) 

(70) 


( 71 ) 


24 



From Equations (69) through (71) and (50) through (5?) we see that 


the spatial statistics associated with particle and rough surface 
scattering are identical. The temporal statistics, however, can be 
considerably different. For rough surface scattering the temporal 
statistics are determined entirely by the source, while for particle 
scattering the temporal statistics are strongly influenced by particle 
motion. 

C can be simplified if we assume the particle velocity is a Gaussian 
h 

random variable. The Gaussian assumption is valid for Brownian motion and 
for turbulent flow which may be found in clouds and smoke plumes. Using 
(65), we have 


fT 


C £ (t) = <exp 


ik • VqT + ik 


v(0 


d£ Av(0 > 


(72) 


where is the average particle velocity an. 4 Av is the fluctuating 

velocity component. The first term in (72) represents the Doppler 

frequency shift associated with the sc cered field. Because we are 

using a direct detection receiver, the Doppler term will have ro effect 

on the received signal and can be dropped. Thus, C may be written in the 

h 

form 


_ x 

C = <e > 
E 


x 


ik 


T 

d£ AvU) 

Jo 


(73) 


Since x is a zero mean Gaussian random variable, we have 


* <r > ■ exp 


H <x2> ) 


2 . 2 

<>: > = k 


n 


(74) 


d4. 


d 2 C Av (c, r' , 2 ) 
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From Equations (69) through (71) and (50) through (52) we see that 
the spatial statistics associated with particle and rough surface 
scattering are identical. The temporal statistics, however, can be 
considerably different. For rough surface scattering the temporal 
statistics are determined entirely by the source, while for particle 
scattering the temporal statistics are strongly influenced by particle 
motion. 

C_ can be simplified if we assume the particle velocity is a Gaussian 
random variable. The Gaussian assumption is valid for Browniar motion and 
for turbulent flow which may be found in clouds and smoke plumes. Using 
(65) , we have 




T " 

= <exp 

ik • VqT + ik • 

d£ Av(£) 



0 J 


(72) 


where v^ is the average particle velocity and Av is the fluctuating 
velocity component. The first term in (72) represents the Doppler 
frequency shift associated with the scattered field. Because we are 
using a direct detection receiver, the Doppler term will have no effect 
on the received signal and can be dropped. Thus, C £ may be written in the 
form 


_ x 

C E = <e > 


x - ik 


T 

d e Av(0 

Jo 


(73) 


Since x is a zero mean Gaussian random variable, we have 


C £ = <e > = exp[~y 


("2 <X2> ) 


2 , 2 
<x > = k 


r~ 


dc 


(74) 
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where C is the velocity correlation function. If the velocity fluctua- 

2 

tions are stationary, <x > can be written as 


2 

< x > 


ni 2 

-2k t 


(1 - e/x) c Av (o 


(75) 


It is interesting to evaluate (75) for the limiting cases where x is large 
compared to the velocity correlation time and where t is small compared to 
the velocity correlation time. If we let denote the velocity correlation 
time, (75) becomes 

2 2 2 

k t <Av > t << T 


<x > 


- 


r 




C iv (t) 


T >> I 


(76) 


Notice that the correlation function C takes on the Gaussian form 

L 

whenever t << x . The velocity correlation time is negligible for 

Brownian motion and in this case C takes on the exponential form which 

h 

gives rise to the familiar Lorentzian spectrum. For particles suspended 

in a turbulent flow, the velocity correlation time >s usually quite long 

( ms) compared to the observation time so that is Gaussian. 

Normally the rotational motion of the particles is negligible compared 

to the translational motion. Either the source coherence function C„ or C_ 

1 E 

dominates the temporal coherence of the scattered field. Consequently, if 

we drop C from (71) and use the Gaussian distribution in (53) for A , 

J finally simplifies to 
A 


a* 



(r,t ) 


l.C (t) C ( t ) exp 

11 L 


1 1 + Y 1 I 2 

2 



(77) 
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The analysis in this section must be modified slightly if the scattering 
mechanism is resonant absorption and reradiation by molecules. For 
resonant scattering, only the incident energy lying within the molecular 
absorption band will be scattered. Therefore, the temporal coherence of 
the scattered field will depend only on the motion of the molecule. The 
effect can be included in an analysis by simply replacing in (77) 

by Ij, the incident intensity lying in the absorption band. If the 
total bandwidth of the transmitted signal is narrower than the absorption 
band, then "f will be equal to 1^. 
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VI. POLARIZATION EFFECTS 


The analysis in the previous sections considered only one polarization 
component of the scattered field. In many cases the scattering process 
will partially depolarize the scattered field so that the total field in 
the observation plane should be written 


A(r,z,t) = A (r,z,t) x + A (r,z,t) jr 
— — x — — y - 


(78) 


where x and y are the unit vectors in the x and y directions. The total 
intensity is given by 


A 


A 


= I (r,z,t) + I (r,z,t) 
x y 


(79) 


Each polarization component contributes a speckle intensity pattern. 

Depending on the characteristics of the scattering medium, there can be 

an arbitrary degree of correlation between I and I . Goodman"* has shown 

x y 

that by using a suitable coordinate rotation, the partially polarized 
field can be decomposed into two uncorrelated polarization components 
of different intensities. If we denote these two polarization components 
by 1^ and 1,, their average intensities are given by 


<1 1> = \ <I T >(1 + f 1 ) 
<1 2 > = j <I T >(1 - P) 


(80) 


where <1 > is the total average intensity and P ij the degree of polarization. 
Because 1| and 1^ are uncorrelated, the intensity covariance function 
can be written as 


C, = Cj (rj ,r 0 , t j , t.,) + Cj (r^r., 




) 


( 81 ) 



The spatial and temporal correlation lengths of 1^ and are identical 

and may be calculated using the procedures outlined in Sections IV and V. 

Thus C and C have the same functional form but differ in magnitude. 

1 2 

This implies 



(82) 


Since the analysis in the previous sections was based on completely 
polarized fields, it is convenient to write (81) in terms of the 
intensity correlation function for complete polarization, i.e., P = 1. 
If we let Cj denote the correlation function for complete polarization, 
then 



(83) 


Using the results in (80) through (83), we find 

C I (£i»r 2 » t i’ t 2 ) = \ C 1 + p2 ) * ( g M 

When the field is completely depolarized, P = 0 in (84) and the intensity 
variance is decreased by a factor of two. This is expected because in 
this case the receiver is actually adding two independent speckle patterns. 

The analysis in Section III can now be used lo calculate the statistics 
of the receiver output by simply replacing the intensity covariance 
function by and the mean intensity by < Ij > * The final results show that 
the mean value of the detected signal is increased to account for the energy 
in the additional polarization component and the speckle noise power is mul- 
tiplied by | (1 + P ) to account for averaging of the speckle patterns in the 
two components. 
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VII. ATMOSPHERIC TURBULENCE 


Atmospheric turbulence affects a propagating laser beam by introducing 

random amplitude and phase fluctuations. This results in a reduction of the 

spatial and temporal coherence of the beam. Very little work has been 

9 10 

reported on the interaction of turbulence and speckle. ’ However, some 
interesting results were derived in the paper by Lee et al.^ which can be 
applied to the satellite-based lidar problem. 

Lee et al. assumed that phase perturbation of the laser beam is the 
dominant effect due to atmospheric turbulence. In addition they assumed 
that the field statistics at the receiver are jointly Gaussian. Utilizing 
these assumptions, the intensity covariance function of the field in the 
observation plane was derived for collimated and focused beams. The 
results are 


Cj (p) 



( focused) 


and 

Cj(p) 

where 



is the intensity variance and 


(coll imated) 


(85) 


( 86 ) 


2 2 -3/5 

= (0.546C zk ) ' . (87) 

0 n 

2 

Py is the phase correlation length and C is the turbulence structure 
parameter. 

The results in Eqs. (85) through (87) apply to the case where the 
initial beam profile at the transmitter is Gaussian and the transmitter, 
receiver and target are all immersed in homogeneous turbulence. 
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Although this does not apply to the geometry of the satellite-based lidar 

where only the target may be immersed in turbulence, the results can give 

some insight into the satellite problem. 

For no turbulence C = 0 and p_ = °° so that (85) and (86) reduce, 
n U 

as expected, to the cases considered in the previous sections. For strong 
turbulence, the intensity correlation length is given by either or R T 
for the focused case and by either p n , R or z for the collimated case. 
In strong turbulence p^ could be small enough so that it determines the 
speckle scale size. 

Physically, the effect can be explained by noting that turbulence 
reduces the beam coherence as it propagates to the target. Consequently, 
the beam spreads out, illuminating an area larger than the diffraction 
limited spot size. As the reflected field propagates back to the receiver 
it is again perturbed by the turbulence, further increasing the effective 
size of the illuminated target. The e_nd result is a reduction in the 
effective speckle scale size which is inversely proportional to the size 
of the illuminated target. 

Similar effects will occur over the satellite-earth propagation path 

Pq will be different from the value given in (87) because the turbulence 

2 

structure parameter is a function of altitude and the turbulence is 

concentrated near the target. For high-alt itude targets such as clouds , 

2 

turbulence effects may be minimal since is significant only near the 
ground . 

Turbulence will also affect the temporal statistics by reducing the 
temporal coherence of the received signal. Because turbulence coherence 
times are relatively long (^ ms), these effects are probably negligible 
compared to the effects of particle motion and transmitter pulse length. 
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VIII. SUMMARY 


In the previous sections the statistics of the signal and noise at the 
output of a satellite-based lidar receiver were evaluated. These statistics 
are summarized in Tables I and II for the case where all the system functions 
are Gaussian. The mean signal, E(s), may be regarded as the average number 
of signal photons detected in a time interval of width 1/B when the receiver 
FOV is equal to the transmitter divergence angle 0^. Similarly, E(n^„) is 
the average number of background noise photons detected in a time interval 
of width 1/B when the receiver FOV is equal to the diffraction limited 
value of \/R. The signal variance, Var(s), and background noise variance, 
Var(n), consist of a shot noise component and a speckle noise. The speckle 
noise power is a function of the lidar system parameters and the characteris- 
tics of the scattering medium. 

A useful parameter for measuring the relative quality of lidar data 
is the signal-to-noise power ratio (SNR) which is defined as 


SNR 


E(s) 2 

Var(s) + Var(n) 


( 88 ) 


Under strong signal conditions the signal speckle noise dominates and the 
SNR becomes 


SNR = 


2M s M T 

(i + r 2 ) 


(89) 


It is instructive to evaluate (89) for some typical lidar parameters, 
As an example we will use the system specification for the proposed cloud 
climatology experiment (Advanced Applications Flight Experiment, GSFC). 
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TABLE I. 


E [ s (t ) + n(t) ] - 


Var [s(t) + n(t)] ■ Mil 

/2 


E(s) + E(n dl ) 



+ j (1 + P 2 ) 


\ / 
Ms)^. 

M A 


+ 


E(n d t > 

Si 



s(t) - 

n(t) - 

n,„ _ 
cU 

R - 

0 T ” 
\ - 

P - 

M - 

s 

«r- 

N - 

S 


signal voltage at receiver output 
background noise voltage at receiver output 
equivalent noise voltage for diffraction limited FOV 
receiver aperture radius 

transmitter divergence angle and receiver FOV 
optical wavelength 

degree of signal polarization, 0 < P < 1 

effective number of spatial correlation cells seen by 

receiver aperture (signal) 

effective number of temporal correlation cells seen by 
receiver electrical filter (signal) 

effective number of spatial correlation cells seen by 
receiver aperture (background noise) 

effective number of temporal correlation cells seen by 
receiver electrical filter (background noise) 
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TABLE II. 


Scattering Mechanism 

M 

s 


2 2 

ire; 

Rough Surface 

l + v 

Molecules, Resonant 

• 

Molecules, Rayleigh 

• 

Particulates a 
Ae roso 1 s 

• 


*r 


4v 


2 


2 2 2 2 2 2 
B B IK + A^/ep 


1 + ^ + 


4v‘ 


'll 

2 


22 22 22 22 
B t B t B (R^ + A Z /9n 
c P r _j 


(v + ; d)‘ 


1 + 


„2„2 
«\ o_ 
T 


b: 2 

1 + -I + V 

2 22 22 
B B^(R Z + X^/Qt) 


1 

12 


i - source coherence time 
sc 

T d ~ *^ nverse Doppler linewidth of scattered signal 

T - transmitted pulse width 
P 

B - receiver electrical filter bandwidth 
B j - receiver interference filter bandwidth 
v - transverse velocity c 5 spacecraft 


go 

§ 
o 


jd r* 

ss 

Eg 

3 * 








Laser 


Receiver 


X - 1.06 pm, 0.532 pm R - 9 cm 

T p - 5 ns (FWHM) B - 100 MHz 

T - 00 V - 10 km/sec 

sc 

0 - 0.5 mrad 

(Half-angle 10 dB 
points) 


Using the above values in (89) we obtain 


3.6 x 10 X ■ 1.06 pm 


bNR = 


1.45 x 10 X * 0.532 pm 


These relatively high SNRs are a consequence of large transmitter 
divergence angle. They illustrate that speckle noise can be maintained 
at acceptable levels by careful choice of the system parameters. 

Although the equations in Tables I and II were derived for the specific 
case where all system functions are Gaussian, the general expression for the 
signal variance can be written in a similar form 


Var (s) 


6E( s) 


+ a t n 

2m s m t 


E 2 (s) 


f90) 


where 

R, (0) 

6 (91) 

/ dTh(f) 

—OO 


j. f/i-iviM 2 

\ 9 ? 

!|J S (°) 1/ d xW(r)p 


(92) 
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In equations (90) through (93) we have neglected the minor effects caused 
by the satellite motion and have used the fact that the mutual Intensity 
function is separable 


(£ » t ) « J s (i)J t (t) 


(94) 


When the receiving aperture diameter is small compared to the spatial 
coherence length, M^, is approximately one. When the receiving aperture 
diameter is large compared to P^, M g can be approximated by the formula 


V 0) 


s [/drW(r) ] 


/ ,2 t^Cr)! 2 

/d r EF 


(93) 


In this case M g is just the ratio of the receiver aperture area to the 
effective area of the spatial component of the mutual intensity function. 

Mj. is also approximately one when the bandwidth of h is small compared to 
the temporal coherence bandwidth of the received signal. When the bandwidth 
of h is large compared to T ^ , M^. can be approximated by the formula 


1 - 


\(3) £ di|3f T (T )|V(T) 


^ ifd-rMi)] 2 |X(0) j 2R p(0) 


(95) 


The effective temporal coherence function of the signal is J^R^. Thus, 

M,p is just the ratio of the effective area of the receiver impulse response 
to the effective area of the temporal coherence function. 

As an example, consider the more realistic situation of a pho»on 
counting receiver where the aperture is annular with an obscuration ratio 
of > and the counting interval is T seconds in duration. Then if f , .1^ 


and P are Gaussian, we have 


m * g 2 a . --T 2 > 

9 2np 2 

c 


2 (1 - 


2 R 2 



( 97 ) 


( 98 ) 


Exact expressions for M g and M^. have been derived for the more 
realistic cases such as annular receiving apertures and higher order laser 
modes. These results will be Included in a £ bsequent report. 
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